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SUMMARY 


The results of some experimental and theoretical studies of the 
interaction of oblique shock waves with laminar boundary layers are 
presented. Detailed measurements of pressure distribution, shear dis- 
tribution, and velocity profiles were made during the interaction of 
oblique shock waves with laminar boundary layers on a flat plate. From 
these measurements a model was derived to predict the pressure levels 
characteristic of separation and the length of the separated region. 


INTRODUCTION 


The interaction of shock waves with boundary layers has been the 
subject of many recent investigations, both theoretical and experimental. 
The theoretical predictions of the whole interaction pattern have, in 
general, been successful only for very weak shocks. The more practical 
problem of strong shocks has remained largely experimental. 

No attempt is made herein to review the many analyses of weak 
interactions. (See, e.g., refs. 1 and 2.) These analyses generally 
involve linearized equations of motion, with or without viscous terms, 
for a flow slightly perturbed from the Blasius flow. 

Another semianalytical approach has been to assume that the boundary 
layer follows a certain pattern and that certain overall parameters are 
sufficient to characterize it throughout the interaction region. To 
this class belong the Pohlhausen procedure of Lees (ref. 3) and the 
integral treatment of Crocco and Lees (ref. 4). A basic assumption of 
this approach is that there is no pressure gradient across the boundary 
layer. This is a reasonable condition everywhere except near the shock 
impingement point and a very useful one because it relates a part of the 
interaction problem to general boundary-layer theory. (This assumption 
is used throughout the analyses in the present work.) However, as Gadd 
points out (ref. 5)> any specific velocity profiles chosen may seriously 
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limit the validity of the analysis in separated flow, and detailed 
experiments are needed to guide the selection of the analytical model. 

The theoretical problem has been eased somewhat by the discovery 
that the flow parameters near separation depend only on local condi- 
tions when the separated region induced by the shock wave is long enough, 
as has been noted in several researches (e.g., refs. 5 and 6 ). This 
forms a basis for the analysis of the present measurements and for the 
derivation of an analytical model. 

The experimental part of this work consisted of measurements of 
pressure and shear distributions and velocity profiles during shock- 
wave — boundary-layer interaction on a flat plate. These measurements 
confirmed that flow parameters near separation are similar. On the 
basis of the observed flow patterns and simple momentum considerations, 
estimates were made of the pressure levels near separation; agreement 
with experiment is good. 

Considerations of conservation of mass and momentum led to a law 
for the variation of length of the separated region with shock strength, 
Mach number, and Reynolds number. The presert measurements verify the 
dependence on shock strength, but more experiments are needed to check 
the Mach number and Reynolds number dependence. 

This investigation was carried out at tie Massachusetts Institute 
of Technology under the sponsorship and with the financial assistance 
of the National Advisory Committee for Aeronautics. 


SYMBOLS 


C 

C P 

C f 

F 

H 


Chapman -Rube sin constant defined by 


mo 


p - p 

pressure coefficient, 2. 



skin-friction coefficient, — 

modified Polhausen-Gruschwitz parameter defined in figure 11 

5 * 

boundary-layer shape factor, — 



height that line u = 0 lies above wall 


Pohlhausen-Gruschvitz parameter 

length of constant-pressure separated region 

Mach number 

pressure 


dynamic pressure. 


* 


^o u o 
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Reynolds number based on distance from leading edge 
Reynolds number based on distance from wall 

recovery factor 
temperature 

component of velocity parallel to plate 
component of velocity normal to plate 

length parallel to plate measured from leading edge of plate 
length from leading edge to end of plateau 

length normal to plate 


r 


i 


specific-heat ratio 

boundary- layer displacement thickness 

height which equivalent undisturbed boundary layer is lifted 
boundary-layer momentum thickness 
dynamic viscosity coefficient 
kinematic viscosity coefficient 
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p density 


t shear stress, \i — 

Subscripts : 

o conditions in undisturbed flow before interaction 

oo conditions at outer "edge" of boundary layer 

B conditions corresponding to Blasius solution 

c conditions which would occur at constant pressure 

d difference between final value and value corresponding to 

incipient separation 

f final conditions after incident and reflected shocks 

pi conditions at plateau 

r reference 

s conditions at separation point 

ST conditions at stagnation point 

w conditions at wall 


EXPERIMENTAL APPARATUS AND TECHNIQUE 


Wind Tunnel 

The experiments were performed in the £- by 8-inch continuous-flow 
supersonic wind tunnel in the Gas Turbine Laboratory of the Massachusetts 
Institute of Technology. This tunnel has a fixed M = 2 nozzle with an 
upward flow direction. A standard schlierer optical system was available 
for observing and photographing the flow. The tunnel can be operated at 
stagnation pressures ranging approximately from 3 to 15 pounds per square 
inch absolute by using a steam-ejector system for evacuation to subatmos- 
pheric levels. However, practical problems in manometer accuracy pre- 
vented successful operation at the lowest available pressure levels. 
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Flat Plate 

As shown in figure 1, a full-span sharp-nosed flat plate was 
installed, approximately on the center line of the test section; details 
of the plate are shown in figure 2. The plate was polished smooth after 
machining in order to avoid any disturbances which could cause transition 
of the laminar boundary layer. Forty-five static-pressure orifices were 
drilled in the plate (using a No. 80 drill) and connected to both mercury 
and Meriam fluid (specific gravity, 2.95) manometers. The orifices to 
measure streamwise distribution of static pressure are spaced one-tenth 
of an inch apart. A thermocouple was installed in the plate surface; 
its temperature was read by means of a standard potentiometer. Stagna- 
tion pressure and temperature were observed similarly, with orifice and 
thermocouple installed in the settling chamber. 


Shock Generator 

The adjustable-angle shock generators developed by Barry, Shapiro, 
and Neumann (ref. 7) for the tunnel were used without changes. 


Velocity-Profile and Skin-Friction Measurement 

The total head in the boundary layer was measured by a small slit- 
mouth tube (fig. 5) which could be traversed streamwise and normal to 
the plate by controls outside the tunnel. The tube height is about 
0.004 inch and the opening is about 0.001 inch. The tube location was 
observed through a cathetometer which was accurate well within 0.001 inch. 
The velocity profiles were computed in the standard manner assuming con- 
stant pressure and total temperature across the boundary layer. This 
method is believed to be sufficiently accurate except in the immediate 
neighborhood of the impinging shock. 

The same tube was used to measure local skin friction by taking 
measurements with the tube resting on the plate. A calibration of the 
tube as a skin-friction meter was made by comparing the tube readings 
with the measurements of an absolute floating-element skin-friction 
meter built specifically for this purpose (fig. 4). The theory for the 
use of small impact pressure probes as skin-friction meters is given in 
reference 8, in which it is predicted that the difference between the 
impact pressure measured by the probe resting on the surface and the 
undisturbed local static pressure is proportional to the 5/3 power of 
the local shear stress, unless the tube is very small or very large. 

The calibration of the tube used in the present work was reported in 
reference 9 and is reproduced in figure 5* It is evident that the 
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5/5-power law is satisfied and no appreciable compressibility effects 
axe observed. 

This calibration was used in the evaluation of the shear-stress 
distributions in the shock-wave — boundary-layer interaction regions 
where pressure gradients exist, although the calibration was made with 
no pressure gradient. (Because of side forces acting on the floating 
element, it furnishes a reliable reference measurement only under zero 
pressure gradient.) Because of the small disturbance region of the 
probe it is believed that the shear- stress measurements are reasonably 
accurate except perhaps where the pressure gradient is very steep. 


Observation of Separation 

The location of separation was taken to be the point where the sur- 
face skin friction, as measured by the total -head tube resting on the 
surface of the plate, is zero. Since the instrument was not able to 
indicate reliable values of negative she?ir, it was often necessary to 
extrapolate the zero position from readings of positive skin friction 
upstream. Although this procedure does not lead to a serious error in 
the location of separation, it can lead uo appreciable error in the 
determination of the pressure at separatr.on because the pressure gradient 
near separation is large. This is probably the major cause of scatter 
in the determination of separation pressure. 

In measurements of the minimum shock strength needed to cause sepa- 
ration, the shock generator was adjusted to provide zero skin friction 
at one point . 


Range of Measurements 

The experiments were performed at several different shock strengths 
up to a static-pressure ratio of 2.4, anc, at different stagnation- 
pressure (i.e., Reynolds number) levels. The shock strengths ranged 
from those small enough not to cause sep«xation to cases involving 
extended upstream influence. The stagnation-pressure range corresponded 
to Reynolds numbers based on distance from leading edge to shock inter- 
action of 1 to 6 x 105. 

Because of the nature of the pressure -level control of the wind 
tunnel, it was not feasible to attempt to take the measurements following 
a rigidly specified parameter -variation scheme. It is believed, however, 
that the results obtained give a reasonat le representation of the flow 
patterns and their characteristics withii. the overall range of Reynolds 
numbers and shock strengths. 
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TYPICAL MEASUREMENTS ON FLAT PLATE 


In figure 6 are shown typical pressure distributions, shear distri- 
butions, and velocity profiles on the flat plate. Figure 6(a) shows 
measurements for a shock not sufficiently strong to cause separation. 

It is seen that the skin-friction coefficient is near the theoretical 
valued upstream of the shock impingement point, falls rapidly to a low 
value near the shock, then again approaches the theoretical value down- 
stream of the shock. Figure 6(b) shows measurements for a case in which 
a separated flow region exists. Again the skin-friction coefficient is 
near the theoretical value both upstream of separation and downstream of 
reattachment. The velocity profiles indicate that in the separated region 
the profile above the reverse flow is essentially the undisturbed constant- 
pressure boundary-layer profile. This will be made a crucial part of the 
subsequent analysis. The pressure distribution shows a constant-pressure 
undisturbed region, a pressure rise to a separation, an additional pres- 
sure rise beyond separation culminating in a constant-pressure separated 
region, a pressure rise up to reattachment, and a further pressure rise 
beyond reattachment culminating in a constant-pressure region corresponding 
to the pressure downstream of the incident and reflected shocks. Fig- 
ure 6(c) also shows measurements for a case involving separation. The 
qualitative picture of the approach to separation and the separated 
region is similar to that shown in figure 6(b). Shortly after reattach- 
ment, however, the skin-friction measurements and the boundary-layer 
profiles show that the boundary layer becomes turbulent. 

Figure 7 shows boundary-layer profiles in the separated region. 

They indicate that the profile above the reverse flow is essentially 
an undisturbed boundary-layer profile. 

Measurements of pressure distributions, shear distributions, and 
velocity profiles, similar to the sample ones given in figure 6 are 
shown in figure 8. 


SIMILARITY OF SEPARATION PROPERTIES 


Figure 9 is a sketch of typical pressure and shear distributions 
which occur when a separated region exists and the boundary layer remains 


^The theoretical skin-friction coefficient is calculated from 
C f = C f-R^ where cp is the Blasius value and C is the Chapman- 


C B* ' X B 

Rubesin constant defined by 

experiments C ~ 0.92. 


M P 

— = C — . For the conditions of these 


Mtdo 


T 
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laminar throughout the interaction region. It is the physical model for 
the following analysis and defines the quantities which will he used. 

Chapman, Kuehn, and Larson (ref. 6) *nd Gadd (refs. 5 and 10) have 
pointed out that, at least when there is ?n appreciable length of sepa- 
rated flow, conditions near separation are independent of the agency 
which induced separation and the separation phenomenon is essentially 
similar; all significant parameters are then functions of Mach number 
and Reynolds number. Thus, the analysis of the shock-wave — boundary- 
layer interaction can be broken up into two parts: (l) The determina- 

tion of similar separation properties, wh:.ch are functions only of 
local conditions, and (2) the relationship between the length of sepa- 
rated region and the shock strength. 

The important parameters involved in separation are the pressure 
rise to the separation point and the pressure rise to the constant- 
pressure separated region (henceforth called the "plateau” pressure 
rise). The minimum pressure rise to cause separation also will be 
analyzed from the point of view of a similar interaction. Chapman, 
Kuehn, and Larson have given an order of magnitude analysis which sug- 
gests that 


AC 


P 



(i) 


In the analysis of Chapman, Kuehn, and Larson the skin-friction coef- 
ficient at the beginning of the separating region Cj> is taken as a 

characteristic skin-friction coefficient for the region, and ACp is 

a characteristic pressure-rise coefficient for the separating region. 
Thus, equation (l) should hold both at separation and at the plateau. 
Gadd’s quantitative calculation (ref. 10) also leads to equation (l) 
for the separation pressure coefficient ar.d determines the constant of 
proportionality for separation. 


Pressure Rise to Separation 

The same law (eq. (l)) for the pressire rise to separation, with 
an estimate of the constant of proportionality, can be obtained from 
simple momentum considerations, using the observation already made that 
the undisturbed boundary layer appears to be "lifted" from the wall. 
Thus, it is assumed that in the approach to separation the boundary 
layer consists of an outer portion, which differs but little from the 
undisturbed layer, and a sublayer having email momentum. A sketch of 
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a typical profile in this region is shown as follows: 



The sublayer joins the outer profile at a height and, relative to 

the outer flow, the undisturbed boundary layer appears to be lifted a 
height e from the wall. If it is assumed that the sublayer has neg- 
ligible momentum, then the force balance in a region of height € gives 


€ 


dp = 

dx 


w 


( 2 ) 


If the pressure rise is taken to be due to the growth of e and small 
perturbations are assumed, the pressure rise is related to the rate of 
growth of the sublayer by 


, 2 de 

,p 


(3) 


Eliminating e between equations (2) and (3) leads to 


2 dx 



( 4 ) 



1C 


An integration by parts gives 


C-n dx = + 

* 2 


Ct>"U) 


Ctj(ti) dr] d| 


Since up to the point of separation Cp" is small except for a small 
region in which C p is small, the integral on the right-hand side of 
equation (5) is small in comparison with C p 2 . Neglecting this integral 
and evaluating equation ( 4 ) at the separal ion point where Cp v = 0 gives 
the separation pressure rise as 


C Ps “ 2 




where 


(° f 0 s 


is the skin-friction coefficient at e at the separa- 


tion point. 


To evaluate 


( Cf e) s 


the form of the sublayer profile must be 


chosen. Equation ( 2 ) implies that the parabolic profile which must 
exist at the wall adequately describes the profile for a small distance 
away from the wall. Thus, an integration of the wall condition 


dp _ / St\ 
dx W/w 


leads to the profile 


T y dp / cL l O 

u = — y + y^ 

t*w ^ My 
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The parabolic profile and the "undisturbed” profile are joined at 
the matching conditions being 



Is €i + gg/fo 
^ 2^ 


ei2 


(9) 



( 10 ) 


T 



, dp 

T 0 = T V + € X 

dx 


(id 


Actually must be a little less than t q because there is an 

adjustment region before the region of fairly constant pressure gradient 
is reached. This adjustment region is ignored. Solving equations (9), 
(10), and (ll) for e gives 



( 12 ) 


At the separation point t„ = 0 and it is found that i t„. 

Thus, the separation pressure coefficient is 


C 


Ps 



(13) 


In figure 10 are shown the separation pressure rises measured in 
the present work. Within the scatter of the experimental points it is 
seen that equation (13) correlates the data reasonably well. As pre- 
viously pointed out, the large scatter is due primarily to the fact that 
small errors in the location of the separation point lead to large errors 
in the separation pressure because of the large pressure gradient near 
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separation. 

For 

C Ps = 

C fo 

p ' 

just given. 

The 


additional comparison, Gaid’s theory (ref. 10) gives 
in agreement with the re salt of the simple reasoning 
data correlated by Chapman, Kuehn, and Larson show a 


slightly lower value; from these data, at Mq = 2, Cp s 



It is interesting to examine the effect of employing the same basic 
assumptions within the framework of the K£rm£n integral equation of 
boundary-layer momentum. For small disturbances from free-stream con- 
ditions this equation can be written 


c f 




+ 0 2 



(im 


Corresponding to the assumption that the boundary layer is lifted by a 
sublayer of negligible momentum, the momentum thickness can be approxi- 
mated by its constant-pressure value. Writing the displacement thick- 
ness as the constant-pressure value plus a deviation from it and writing 
instead of equation (3) 


, 2 a(A6*) 

,p p dx 


where A6* is the deviation from the constant-pressure displacement 
thickness, gives on substitution into equation (l4) the following 
equation: 


% “ C fv = ( H c + 2 - Mo 2 ) Q c ^ + f ^ f x °P (l6) 

x o 


vhere the subscript c denotes the constant-pressure value. As before, 
an integration by parts is performed and then an evaluation of equa- 
tion (l6) at the separation point, where Cf = 0, gives 
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+ 2 - Mo 




(17) 


Now, if the boundary-layer separation mechanism simply involves 
the action of an adverse pressure gradient, then perhaps an estimate 

( dC p\ 

of \0 C can be obtained from the extended Pohlhausen calculation 

\ ax/s 

( e.g ref. ll). The boundary layer is then described by the Pohlhausen- 
Gruschwitz parameter K, which at separation can be written 


K s 


Pco 9^ ^00 ^.Cp j 
P v C 2v ro dx / s 


-0.157 


(16) 


Equation (l8) leads to the approximate result 



(19) 


or, since 0 C = xc f , 


dC T 


x 


0.713 


dx / s 7-1 2 

1 + “TT 


( 20 ) 


Note that although equations (19) and (20) have the form of the incom- 
press ible equations, they actually are valid for compressible flow 
because of the use of the Chapman-Rube sin constant C. 

In equations ( 19 ) and (20) r is the recovery factor. Equation (20) 
is subject to experimental verification. Data from the NACA, the University 
of Southern California, the Massachusetts Institute of Technology, and the 
National Physical Laboratory for a Mach number range of 2.0 to 2.55 and a 
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Reynolds number range of 76,000 to 1,150,000 in which separation was 
induced by a step, a corner, and a reflecting shock wave are plotted in 
figure 11. Considering the difficulty in measuring an experimental pres- 


sure gradient, these data do indicate that tie product 



at separa- 


tion is constant at a given Mach number but lhat the constant is less than 
half the Pohlhausen-Gruschwitz constant. Adjusting the constant and 
putting equation ( 19 ) into equation (17) give the following equation 
for the separation pressure coefficient: 


C 


Ps 




+ l 




( 21 ) 


From the work of Chapman and Rubesin (ref. 12), the following equation 
can be derived for 7 = 1.4 and Prandtl number = 0.72: 


H 


c 



1 + 0 . 



( 22 ) 


The subscript B denotes the Blasius value. Therefore, for M 0 = 2, 
equation (21) shows that 


C 


Ps 



( 23 ) 


The reasonable agreement between equaticns ( 15 ) and ( 25 ) indicates 
that the integral momentum equation may be satisfactory even in strong 
pressure gradients in supersonic flow if it is applied judiciously. The 
failure of the Pohlhausen-Gruschwitz parameter to describe the pressure 
gradient at separation is due to the inadequate representation of the 
boundary-layer profile; the integral equatior is not at all involved in 
this calculation. 

It is expected that the correction to the Pohlhausen-Gruschwitz 
parameter depends upon the Mach number. Lacking an extension of the 
correction to Mach numbers other than 2, the simpler equation (13) 
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should, be used at other Mach numbers. The agreement between equations ( 13) 
and ( 23 ) will be taken as justification for reasoning from the integral 
momentum equation. 


Pressure Rise to Plateau 

Beyond the separation point the analysis is more difficult. Clearly, 
the parabolic approximation for the inner boundary-layer profile cannot 
hold throughout this region both because it would degenerate to an incor- 
rect linear one when the plateau is reached and because the inertia forces, 
far from being negligible, must balance the shear forces when the plateau 
is reached. However, it is expected that the inertia forces and the 
curvature of the pressure distribution will come into play about the 
same time; since they are of the same order and make contributions of 
opposite sign, an order of magnitude can be estimated by neglecting 
both of them. With this approach 




(24) 


This is then really the same as the condition prior to separation. Now, 
c f e must be less than the maximum value of c^ anywhere in the profile 

at the beginning of the plateau, and the profiles show that this maximum 
is less than, say, 0.9c£ . If a symmetrical backflow is assumed at the 

plateau, the following approximate upper limit is obtained for the 
plateau pressure: 




or 


(25) 


c P pl < 1 -9 c p s (26) 

Another estimate of order of magnitude can be made by adopting the 
point of view that beyond separation the line u = 0 forms a wall around 
which the boundary layer is deflected. Then the integral momentum equa- 
tion (l4) can be applied above u = 0. If the line u = 0 is lifted 
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a height h from the wall, then the following equation analogous to 
equation (l6) is obtained: 


c fv 


-( 


H c + 2 - M ( 


o l® c 


dC 


dx 


P + P 

2 dx 


,x 


Cp dx 


dCj 

dx 


n x 


i 


dh , 

— dx 
dx 


( 27 ) 


Nov the pressure distribution in the vicinity of the separation 
point has the following qualitative form: 



At the beginning of the interaction region, between points 0 and A, 
there is a decrease of wall shear and an increase of pressure gradient 
which occurs with negligible pressure rise ontil the pressure gradient 
approaches the separation value. Between points A and B the pressure 
increases at almost constant pressure gradient. Beyond point B there is 
a decrease in pressure gradient, and the reverse flow momentum goes from 
a negligible value to its equilibrium value at the plateau. 

If most of the pressure rise to the plateau occurs before the 
reverse flow momentum becomes important (i.?., before point B which 
would give a lower limit for the plateau pressure), then 


dCp 

dx 




(28) 


When the plateau is reached the plateau pressure coefficient is due 
primarily to the angle which the line u = 0 makes with the wall, and 
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u> 


the boundary layer above u = 0 must have thinned down between the 
separation point and the plateau so that dS*/dx is again small. It 
is expected that the decrease in pressure gradient between point B and 
the plateau is accompanied by an increase of shear stress at u = 0 in 
the same manner as the increase in pressure gradient between points o 
and A is accompanied by a decrease in wall shear. If it is assumed that 
the deviation of the skin-friction coefficient at u = 0 from the 
constant-pressure value is the same at points A and B, then if the pres- 
sure at point B is identified with that at the plateau and the addi- 
tional assumption is made that c p ^ == -c f^ when the plateau is approached, 

equation (27) leads to 

Cp pl > Cp s ^ ( 2 9 ) 


Taking the mean between the estimates of equations ( 26 ) and (29) as the 
approximate plateau pressure coefficient gives 



1.65Cp 


S 


Using equation (13) to give C pg 


there results 


(30) 



(31) 


The data plotted in figure 10 show good agreement with the results 
of equation (31) • The correlation of data in reference 6 gives, for 


“o = 2, C ppi = 2 



which is about 13 percent lower than the value 


given by equation ( 3 l)* On the other hand, the correlation indicates 
that 1*68 times the correlated Cp g , in agreement with 

equation (30). 


In figure 12 is plotted the variation with Mach number of the 
separation and plateau pressure coefficients as given by equations ( 13) 
and (31) and by the correlation of data of reference 6. Also shown is 
the data of the present work. Considering the scatter of the experimental 
data it can be concluded that the separation and plateau pressure coef- 
ficients are reasonably well predicted by equations (13) and (31). 
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Pressure Rise for Incipient Sept-ration 

For shock strengths just great enough to create zero vail friction 
at one point, the observed, flow pattern is as shown in figure 15* In 
the vicinity of the shock the pressure gradient on the surface is main- 
tained reasonably constant by the external shock configuration, although 
the self- induced gradient would be negative as the boundary-layer edge 
becomes convex. It can be argued, however, that the local effect of the 
shock system cannot extend very far from the impingement point since the 
pressure jump is progressively canceled by the expansions when the sur- 
face is approached. This view seems to be supported by the experiments. 
Thus, the pressure rise is due primarily to the self-induced effect, and 
the pressure rise to separation should be about the same as that for 
separation well upstream of the shock. If symmetry of thickening and 
thinning of the boundary layer is assumed, there results 

( C p) incipient = 2C Ps (32) 

The data shown in figure 10 support rather well this simple point of 
view. 


LENGTH OF SEPARATED REGION 


When the flow is completely separated, the height of the separated 
region depends on the mass of fluid forced intc the backflow by the 
overpressure left available from the main shod after subtracting the sum 
of the pressure rise from undisturbed conditions to the plateau and the 
pressure rise associated with the reattachment process. Now, all the 
fluid which is reversed must come from above the line u = 0 as shown 



Plateau 
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If the boundary-layer equations are used again, the integral of the 
momentum equation is 


pu — + pv — + 
i dx dy 


&W . r * 

d X j J ($y 


dx 


(33) 


This integral is to be evaluated along the line u = 0 from the maxi- 
mum height of the separated region to the reattachment point. At the 

beginning of the reattachment process — must be of order pv , [ — ) 

dy pi \dy/ 

In the approach to reattachment 


- \°y/ p i 

| v j » | v pl| Por l° n S plateaus, and 


the pressure gradient increases rapidly. It is therefore expected that 
Tor a while the opposite effects on — of the adverse pressure gradient 

°y 

and the suction through u = 0 would tend to cancel each other, and hence 

r St 


it seems fair to neglect 


^y 


dx. Thus, 


Ap, « 


r du 

J PV ^y 


dx 


(3U 


where Ap^ is the driving pressure for the reverse flow. Taking averages 
gives 


Ap a ~ _p v v ^ ^ 


(35) 


By continuity, ~P W V A* = P w uh max , and thus equation (35) becomes 


^W 

^Pd ~ = Py^max 


(36) 
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A momentum balance in the separated region at gives approximately 


_2 dh 3 _2 n 

- U to = 2 Si = T u=0 " 


where the reverse flow has been assumed to be symmetrical. 

Eliminating u between equations ( 36 ) ancl (37) and writing 


+ - C. , l 


2 Ppi 


where Z is the length of the plateau gives , for a value of h^ which 
is small in comparison with h max , 

Cf 5 qoPvP c p pl * 2 = tV^c Pc (39) 

Replacing Cp ^ by its value from equation (31) and letting x r be 
any reference length, there results from equation (39) 


/ \ 5/2 

i; - 2 - 1 ( 1 + rM o 2 J P' 1 ^ c f r 2 cy" 5 / 2 c fo - 1 / 4 ACp. 


where c^ is the constant -pres sure skin- friction coefficient at the 
I r 

reference station x r . If the reference length is taken as the 

length from the leading edge to the end of the plateau, then c- = c- 

Ir 
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In the absence of mixing, Cj> ~ . Accepting this value for the 

x 2 

moment, equation (40) leads to 



The important results of equation (42) are that the length of the 
separated region depends strongly on the driving pressure and only 
weakly on the Reynolds number. 

The driving pressure is not too easily determined from experiment, 
and it would be useful if it were not necessary to do so. A plausible 
choice for the driving pressure is the difference between the final 
pressure and the pressure for incipient separation. Thus, 


AC Pd C P f " ( C p) incipient = °Pf " 2C Ps " °Pf " 1,21C Ppl ^ 


The last form is the most useful, since CL, and C n are easy to 

*f v pl 

measure and the other quantities are not. With this choice of ACp d , 

the lengths of the separated regions are plotted in figure l4. The 
linear relationship between the function of length of the separated 


region and the driving pressure is clear 

able in these tests^. However, the constant of proportionality is 

greater than that of equation (42) and the line does not go through 
the origin. Both of these results are not unexpected. For short lengths 
of separated region the reverse flow is not well established, the wall 
shear must be close to zero, and, hence, the effective must be much 


1/4 


is not a major vari- 


less than the local constant-pressure value. Consequently, for small 
lengths of separated region the length should grow more rapidly than 
for larger lengths in which the backflow is well established. The con- 
stant of proportionality is larger than that predicted by equation (42) 
because c P really is somewhat less than c f . Using the experimental 

Xj 

slope and intercept, equation (42) becomes, for the longer separated 
regions. 
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The new constant of proportionality implies that ^ 0.82c^ which 

x 7 

is quite a reasonable figure. Thus, the experiments seem to support 
the proposed model. More experiments are needed, however, to verify 
the Mach number and Reynolds number dependence. 


CONCLUDING REMARKS 


An experimental and theoretical analysis of the interaction between 
a shock wave and a laminar boundary layer has been made. The qualita- 
tive picture of the interaction is as follows: 

The adverse pressure gradient due no a shock wave causes the wall 
shear to decrease upstream of the shock wave. When the pressure rise 
reaches a certain threshold the wall shear is brought to zero. For 
larger pressure rises a separated region occurs. 

The experiments show, as also is reported in other researches, 
that when the separated flow region is sufficiently long the pressure 
levels associated with separation depend only on local conditions. It 
is observed from the boundary-layer profiles that the profiles are 
modified primarily near the wall, the outer flow remaining essentially 
an undisturbed boundary layer. An analysis of the separation pressure 
levels based on this observation gives reasonable agreement with the 
measurements. It is also shown that the integral momentum equation can 
be used to predict the separation pressure levels provided that an 
adequate measure of the pressure gradient at separation can be made. 

An analysis of the length of the separated region shows that it is 
almost proportional to the pressure rise over and above that needed to 
induce incipient separation and that it varies only weakly with the 
Reynolds number; the length is inversely proportional to the eighth 
root of the Reynolds number. The measurements also support this analysis 


Massachusetts Institute of Technology, 
Cambridge, Mass., May 1, 1957* 
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Figure 2.- Details of flat plate 




Figure 3*- Stanton tube configuration 











Figure 6.- Sample pressure and wall shear distributions and velocity 

profiles . 
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Figure 7*- Velocity profiles in separated flow 
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Figure 9 >- Sketch of flow pattern with separated region. 







SKIN-FRICTION COEFFICIENT AT START OF INTERACTION, c f 

'o 

Figure 10.- Separation, plateau, and incipient separation pressure 

coefficients . 










Figure 14.- Length of separation region. 
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